Symmetric geometric measure and dynamics of quantum discord 
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A symiiietric measure of quantum correlation based on the Hilbert-Schmidt distance is presented 
in this paper. For two-qubit states, we simplify considerably the optimization procedure so that 
numerical evaluation can be performed efficiently. Analytical expressions for the quantum correlation 
are attained for some special states. We further investigate the dynamics of quantum correlation 
of the system qubits in the presence of independent dissipative environments. Several nontrivial 
aspects are demonstrated. We find that the quantum correlation can increase even if the system 
state is suffering dissipative noise. Sudden changes occur, even twice, in the time evolution of 
quantum correlation. There is certain correspondence between the evolution of quantum correlation 
in the systems and that in the environments, and the quantum correlation in the systems will be 
transferred into the environments completely and asymptotically. 
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I. INTRODUCTION 

Quantum systems can be correlated in ways inaccessi- 
ble to classical objects. There are quantum states that 
cannot be prepared with the help of local operations and 
classical communication, or that cannot be represented 
as a mixture of product states. These states are called 
entangled states. Entanglement is certainly a kind of 
quantum correlation, and moreover is by far the most 
famous and best studied kind of quantum correlation 
[1]. One reason for this situation is the fact that quan- 
tum entanglement plays an important role in much of 
the research of quantum information science [2]. How- 
ever, quantum entanglement is not the only aspect of 
nonclassicality of correlations and does not account for 
all nonclassical properties of quantum correlations. A 
typical example is that quantum nonlocality can arise 
without entanglement [3] . More importantly, although it 
is well known that entanglement is essential for certain 
kinds of quantum-information tasks like teleportation [1] 
and super-dense coding [5], there is no definite answer 
as to whether all quantum algorithms that outperform 
their best known classical counterparts require entangle- 
ment as a resource. Indeed, there are several instances 
where we see a quantum improvement in the absence or 
near absence of entanglement (see e.g. [(>-!)]). In par- 
ticular, no entanglement is present in the computational 
model referred to as "the power of one qubit" with the 
acronym DQCl [lU, 11]. Despite this, the mixed sepa- 
rable states can create an advantage for computational 
tasks over their classical counterparts [11, 12]. 

On the other hand, many works have been devoted to 
understanding and quantifying the quantum correlations 
beyond entanglement [13-27]. As the total correlation 
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can be split into a classical part and a quantum part 
[1.1], various measures of quantum correlations are pro- 
posed by considering different notions of classicality and 
operational means to quantify quantumness. Amongst 
them, quantum discord [14] has attracted much atten- 
tion. Quantum discord has been defined as the mismatch 
between two quantum analogues of classically equivalent 
expression of the mutual information. It can be also ex- 
pressed as the difference between total correlation, mea- 
sured by quantum mutual information, and the classi- 
cal correlation defined in [13]. The notion of quantum 
discord goes beyond entanglement: separable states can 
have nonzero discord. In particular, it is believed that 
quantum discord is the figure of merit for the DQCl 
model of quantum computation [2s]. Generally speak- 
ing, quantum discord plays an important role in quantum 
information processing [2!i]. 

In order to obtain the classical correlation (and thereby 
the quantum discord) of a bipartite quantum system, one 
has to perform measurement on one subsystem to extract 
the information about the other subsystem, i.e., locally 
accessible information. Hence this sort of one-sided mea- 
surements implies that classical correlation and quantum 
discord are not symmetric under the permutation of sub- 
systems. As a measure of correlation, whether classical 
or nonclassical, one would expect that it is symmetric. 
Some symmetric measures of correlations have been pro- 
posed in Refs. [18, 19, 22, 30]. The key step in obtaining 
these measures is to consider the classical-classical (CC) 
states. The state p is called a CC state, if it can be ex- 
pressed as the mixture of locally distinguishable states, 
namely, Pcc = J2i,jPij\'i'){^\ ® lj)0'L where Pij is a joint 
probability distribution and local states \i) and \j) span 
an orthonormal basis. The set of product states is the 
subset of the set of CC states, while the set of CC states 
in turn is the subset of the set of separable states. A CC 
state is in fact the embedding of a classical probability 
distribution in the formalism of quantum theory and as 
such has no quantumness. Given a state p, one can find 



the closest CC state, pcc^ to it. It is natural to regard 
this minimal distance as the measure of quantum cor- 
relation. Such a measure has a transparent geometrical 
meaning. The distance can be measured with trace dis- 
tance, Hilbert-Schmidt distance, Bures distance, relative 
entropy and so on. In [21'], the authors used the relative 
entropy as the distance measure to provide the unified 
view of quantum and classical correlation. 

We present in this paper a symmetric geometric mea- 
sure of quantum correlation based on Hilbert-Schmidt 
(HS) distance. By applying von Neumann measurement 
on each subsystem, any bipartite state will become a CC 
state, which we call a measurement-induced CC (MICC) 
state. Given a bipartite state p, we define the quantum 
correlation, denoted by G(p), as the squared HS distance 
between p and the closest MICC state. The evaluation 
of G{p) requires an optimization procedure over the set 
of all local von Neumann measurements on two subsys- 
tems, and thus attacking the general case is a formidable 
task. However, for 2-qubit states, we are able to sim- 
plify the optimization over two-sided measurements to 
that over one-sided measurements. In other words, the 
number of measurement parameters over which the opti- 
mization procedure is performed is reduced from four to 
two. Thus we are able to evaluate the geometric measure 
of quantum correlation by the efficient numerical method. 
We further find the exact analytical expressions for some 
special states. 

Moveover, we will discuss the dynamics of quantum 
correlation, quantified by the geometric measure G, 
when quantum states undergo a noisy channel. It has 
been shown that the evolution of quantum correlation, 
whether in terms of quantum entanglement or of quan- 
tum discord, may behave in a "sudden" way. Quantum 
entanglement can evolve to sudden death or birth [31]. 
The decoherence regime of quantum discord may change 
suddenly to that of classical correlation [■)2]. Also, the 
dynamics of quantum discord has attracted much atten- 
tion [33-35]. This "sudden" behavior has not been thor- 
oughly understood yet. It is natural to enquire as to 
whether and how such a geometric measure of quantum 
correlation will exhibit some sort of sudden changes. To 
this end, we consider the case that a 2-qubit state is ef- 
fected by the action of two independent non-unital chan- 
nels, e.g., amplitude damping (AD) channel, and calcu- 
late analytically G(t), the time evolution of the geometric 
measure developed in this paper. We observe the non- 
trivial phenomenon that the function G{t) may be neither 
monotone nor smooth. In other words, during the time 
evolution, the quantum correlation can increase under 
the influence of AD channel, and the rate of evolution 
can exhibit sudden changes at some critical times. De- 
spite these novel aspects, the asymptotical behavior of 
G{t) is in accordance with the coherence decay. In addi- 
tion, we see that the quantum correlation in the system 
qubits is completely transferred to the environments after 
sufficiently long time. 



II. DEFINITIONS AND NOTATIONS 

This section is a prelude providing the definitions and 
notations that will be used throughout the whole paper. 



A. Hilbert-Schmidt distance 

For any 2-qubit state p shared by two parties A and 
B, let's define a 4 x 4 matrix R as R — (i?pjy) with the 
elements R^i, given by R^^, — Tr[p(crp ® a^)] for p,,v = 
0, 1, 2, 3, where ctq is 2 x 2 identity matrix and ai {i — 
1, 2, 3) are usual Pauli matrices. We write R as 



R = 



1 y 



(1) 



where x and y are the Bloch vectors (in row form) of 
reduced state p"^ and p^ respectively, and T is a 3 x 3 
matrix and usually called correlation matrix. The super- 
script T denotes transposition. 

Hilbert-Schmidt (H-S) norm (also called Frobenius 
norm) is defined as ||M||hs — \/Tt{MM^) for any 
bounded operator M on complex-dimensional Hilbert 
space. It follows that the H-S distance between two op- 
erators M and N is given by 



D 



HS 



(M, N) = yjTT[{M - N){M'< - N't)]. (2) 



In the following, we will omit the subscript HS. 

Let Rp and Rr be the R matrix associated with 2-qubit 
states p and r respectively. The H-S distance between p 
and r can be expressed in terms of the elements of R 
matrix: 



D'^iP,T) = -[\xp-x 



■r\^ + \yp-yr\' + \\Tp-Tr\\']. (3) 



Now suppose that we perform local von Neumann mea- 
surements on both qubits A and B. The measurement 
operators for A and B are given by 

respectively, where k and £ are unit vectors in three- 
dimensional real space. After measurements, we obtain 
a MICC state x, that is, 



X 



J2 (Hf®Hf)p(Hf®H 



(4) 



i,j=+,- 



Let K — k'^ k and L = l'^ £. Both of them are 3x3 real 
symmetric matrices. The R matrix of x can be written 
as 



^^~ ^Kx^ KTL 



(5) 



From (3), the squared distance between p and x is given 

by 

= ^ [\x-xK\^ + \v- yL\^ + \\T- KTlf] 



\>c 



= 4i'^' + y' 



|r||2 - [xKx'^ + yLy'^ + Tr 



{TLT^K)]}, 



where x'^ = x- x and so forth, and in the last hne we have 
used K'^ — K and L^ = L. Let X = x"^ x and Y = y'^ y. 
We rewrite D^{p,x) as 



i^2(p,;^)=l[Tr(x + r + rr^) 



- Tr (XK + YL + TLT^K) 



(6) 



If only one qubit, say A, is measured, the resulting 
state is a classical-quantum (CQ) state. A general CQ 
state has the form of X]j-PiN)(*l ® Pi- I^ the case we 
are considering, the measurement-induced CQ state is 
expressed as 

p-^ = {U^ (E)l)p (U^ (g)l) + {U^ (g)l)p (U^ (E) 1). (7) 

The corresponding R matrix is 



Ro 



1 y 

Kx^ KT 



(8) 



It follows that the squared distance between p and p~^ is 
given by 

D^{p,p-^) = ^[TT{X + T'^T)-Tr{XK + TT'^K)]. (9) 

If only qubit B is measured, the squared distance 
D^(p,p^) can be defined similarly. 



B. Quantum channel 

We will give only a brief introduction to the quantum 
channel, and refer the reader to books, say, [2, 'Ui] for a 
detailed analysis. 

A quantum channel is a trace preserving completely 
positive (TPCP) map. A linear map £ is completely 
positive if and only if it is of the form 



^' = J2K^pKJ, 



(10) 



where Ki are the operators on the state space of the sys- 
tem and are known as Kraus operator. Eq. (10) is called 
operator-sum representation of a quantum operation. A 
completely positive map is called trace-preserving if and 
only if 



J2k}k,^i. 



(11) 



If a quantum channel leaves the maximal mixed state 
invariant, it is called a bistochastic map or unital channel, 
that is, ^^KiKj — 1. Otherwise, the channel is non- 
unital. 

The action of a quantum channel on a quantum system 
can be described by a unitary evolution on an extended 
system, the system S plus the environment E, followed 
by a partial trace over E. This description is called the 
unitary representation of the quantum channel. This rep- 
resentation is not unique since many different unitary 
evolutions will lead to the same effect of the channel. 

Now let us consider a concrete quantum channel, am- 
plitude damping (AD) channel. The AD channel is used 
to describe the evolution of the system state in the pres- 
ence of a dissipative environment. In this process the sys- 
tem interacts with a thermal bath at zero temperature. 
This process could be described as spontaneous emission 
of a two-state atom (system S) coupled with the vacuum 
modes of the ambient electromagnetic field (environment 
E) which leads the atom state to the ground state (see, 
e.g., [37] for more in-depth discussion). 

Considering the behavior of a two-level atom in a A^- 
mode cavity, the interaction between them gives rise to 
the phenomenological map 

|0^)|0^) ^ |G^)|0^), (12) 

|1^)|0^) -^ 7{t)\l')\0'') + Vl-7^(i)|0^)|l^). (13) 

Here, 10"^) is the ground state and jl'^) is the excited 
state of the atom, while |0^) is the vacuum state of the 
cavity and |1^) describe the cavity state with only one 
excitation distributed over all modes. The amplitude ^{t) 
converges to j{t) — exp(— Ki/2) in the limit of A^ — > 
oo. The parameter k is usually called coupling strength. 
Then the atom and environment evolve as an effective 
2-qubit system. 

From (12) and (13), we can write the Kraus operators 
of the AD channel: 



Kn 



1 
7 



Ki 



yi 





r 







(14) 



The evolution of the system state is then given by p^ (t) = 

Let us further consider the case of a 2-qubit system 
being affected by their independent dissipative environ- 
ments. We assume that at time i = the system-plus- 
environment state is described by 



ABA'B' 



P 



„AB 



^A'nB 



0^0") (0^0 



B'\ 



(15) 



where |0 0^ ) is the vacuum state of two environments. 
By means of Eqs. (12) and (13), we can work out 
pABA'B'f^^^^ the total state of ABA'B' at time t. Trac- 
ing out environments A' and B' , we obtain the system 
state p^^ (i) , which can be expressed in the operator-sum 
representation. 



„AB 






^K,)p^^{kI^k]). (16) 



On the other hand, the reduced environment state is 
given by 



A'B'/,\ rp ABA'B' i,\ 



(17) 



III. 



GEOMETRIC MEASURES OF QUANTUM 
CORRELATION 



As stated in the Introduction, we can define the quan- 
tum correlation in a 2-qubit state p as the geometric dis- 
tance, measured by HS norm, between p and the closest 
Mice state x- We will present below a detailed analysis 
in this direction. 



A. One-sided measure 

To begin with, we discuss a simpler case in which only 
one-sided measurement is performed on qubit A. The 
geometric measure of quantum correlation based on such 
a sort of one-sided measurements has been proposed in 
[23] and discussed in [38] . We recover these results below 
for completeness. 

Given a 2-qubit state p, one-sided von Neumann mea- 
surement on A induces the CQ state p~^ given by (7). 
With the squared distance D^{p,p'^) given by (9), the 
geometric measure of quantum correlation is defined as 

G^(p)=minZ?2(p,p^), (18) 

where the minimization is performed over all von Neu- 
mann measurements on qubit A. 

For convenience, we use Dirac bra-ket notation to ex- 
press the vectors in real three-dimensional space, that is. 



v'^ — \v) 



and V — 



We will use alternatively both 



notations in the present work. 

It can be seen from (9) that to find the minimal value 
of D^{p, p~^) is equivalent to calculate the maximal value 
oiTT:{XK+TT'^K). By noting that Tr(Xi^-Hrr'^if) = 
{k\{X + TT'^)\k), we have that among aU unit |fc) G E^, 

(fc|X-|-rr^|fc)^ax = A,„ax, (19) 

where Amax is the largest eigenvalue oi X + TT^ . Thus 
we have 



G^{p) - -[Tr(X + TT'^) - A^ax]. 



(20) 



It is exactly the result given in [23] . Similarly, if qubit B 
is measured, we have 



G-(p) = ^[Tr(y + r^T)-A.„ 



(21) 



with Kniax the largest eigenvalue of y + T^T. Generally, 

In [_ .], the authors have used general, rather than 
measurement-induced, CQ state in the derivation of (20). 
Subsequently, it is pointed out in [38] that such a geomet- 
ric measure is essentially measurement-oriented, meaning 
that the optimal CQ state in the most general sense is 
indeed the sort of measurement-induced CQ state. 



B. Two-sided measure 

Now we proceed to the case of two-sided measure- 
ments, which will lead us to a symmetric measure of 
quantum correlation. Two-sided measurements will re- 
sult in Mice states. The squared HS distance D^{p,x) 
between the state p and the corresponding MlCC state 
X has been given by (6). We define the geometric mea- 
sure of quantum correlation, based on two-sided measure- 
ment, as 



G{p) 



min D^{p,x), 



(22) 



where the minimization is performed over all x, or equiv- 
alently, over all measurement directions {k,£}. As a CC 
state does not incorporate any quantum correlation, the 
two-sided quantifier G can be regarded as a more strict 
geometric measure of quantum correlation. 

By referring to (6), we need only focus on the maximal 
value of Tt{XK + YL + TLT^K) for all allowed K and 
L. Let's define the following two matrices: 



M = X + TLT^ + (£|r|£)l3, 
N = Y + T'^KT + {k\X\k)l3, 



(23) 
(24) 



where I3 denotes the 3x3 identity matrix. It follows 
that 

Tr{XK + YL + TLT'^K) = {k\M\k) ^ {(.\N\e). (25) 

Then the problem is to calculate the maximal value of 
{k\M\k) or {£\N\i) over all unit vectors |fc), \t) € E^. 

By noting that the eigenvalue of M is equal to the sum 
oi {e\Y\£) and the eigenvalue oi X + TLT'^, we see that 
the eigenvalue of M is the function of £, which we denote 
by Am(^) — Am(^i,^2,^3)- Let A™/'^ be the maximal one 
over all unit vector £. It follows that the maximal value of 
{k\M\k) is exactly the A^^''. Similarly, let AAr(fc) be the 
eigenvalue of N and A™'^'^ the maximal one. The maximal 
value of {i\N\i) is \'§^''. Also, we have A^^f'^ = A^'^'^. The 
directions of the optimal measurement, denoted by k °p* 
and £ °p* , are the eigenvectors of M and N referring to 
the eigenvalues A"'*'' and A^^'^, respectively. Thus we 
see that the problem of two-sided optimization is indeed 
reduced to that of the one-sided optimization. 

To proceed further, we will use the following lemma, 
which can be proved by direct calculation. 

Lemma 1. For any two vectors \a) and \b) (not nec- 
essarily normalized) in M.^ , the largest eigenvalue of the 
matrix \a){a\ + \b){b\ is given by 

l[a^ + b^ + ^{a^-b^)^ + 4{a\by], (26) 



A 

with a} — (a I a) and b^ 
malized eigenvector is 

1 



{b\b) . The corresponding nor- 



|A) 



(a^ -b^ + ^(a2-62)2 + 4(a|6)2) 

2b 



a\b) \eb) 



where jco) and |ef,) are the unit vector along \a) and \b) 
respectively, and Nx is the normalization factor. 

Let \k') = T^\k) and |f ) = T\i). It follows from the 
Lemma 1 that 



Am(0 



XN{k] 



2{£\Y\£)+x'^ +e'^ 
2{k\X\k)+y^ + k'^ 



(27) 



(28) 



The rem.aining problem, is to find the maximal values 
or XJ^^^. The geometric measure is then given by 



A 



M 



G{p) = 



1 



y 



\Tr-\ 



M(N) 



(29) 



Referring to (27), we see that Xm{£) involves two param- 
eters, that is, two angles indicating the direction of £. It 
is not difficult to attain the A™/'^ by efficient numerical 
method. For some special states, the exact results are 
available, which will be stated as follows. 



1. States with two-sided maximally mixed marginals 

In this case, x = y = 0. We can always trans- 
form the correlation matrix T into the diagonal form 
diag(ii,t2,i3) by local unitary operations. Eq. (27) be- 
comes Xm{£) = f^ = X]i=i ^f-^f ■ The maximum is given 
hyX'^r = max{tltltl}. 

In this case there is no difference between the one-sided 
measure G^{p) and the two-sided measure G(p). 



The exact expression of G(p) for any X state is in fact 
available. However, the derivation is too lengthy and too 
tedious. So we present the results referring to a restricted 
class of X states, that is, the X states with 1^3 1 = |2/3|,i.e., 
with the identical local purity. The concrete analytical 
results are presented in Appendix. 



IV. DYNAMICS OF QUANTUM 
CORRELATION 



Now we are in the position to analyze the dynamics 
of quantum correlation, in terms of the geometric mea- 
sure developed in this paper. We focus on a non-unital 
channel, AD channel. It is assumed that the identical 
AD channel is applied on each qubit independently and 
simultaneously. In order for the analysis to be precise, 
we take the system's initial state p"^^ as the X state with 
the same local purity. The time evolution of p is given 
by (16). The environment state p^ ^ (i) can be attained 
by considering the total state p"^^^ ^ (i) and then trac- 
ing out A' and B' . We see that at any time both p^^{t) 
and p"^ ^ {t) are the sort of the X state with the identical 
local purity, which allows for an analytical computation 
of the geometric quantum correlation G. 

Let us consider the following two examples. 

Suppose that two 2-qubit X states p^^ and p"\^ are 
given by, in terms of the elements in the R matrix. 



pr 



and 



Pi 



X3^y3 = 0.7949, 

T = diag(0.4705, -0.5277, 0.8947), 



3:^3 = 2/3 - 0.6479, 

T = diag(0.3926, -0.0772, 0.0360), 



(30) 



(31) 



2. States with one-sided maxim,ally mixed margma 



In this case, only one reduced state, say p , is max- 



imally mixed. Then a; = 0. It follows from (27) that 
Xm{£) = {£\TT^ + Y\£) and the A^^'^ is given by the 
largest eigenvalue of the matrix TT'^ -I- Y . On the other 
hand, if p^ is maximally mixed, we can consider X^ik) 
with y — 0. Similar results are easily obtained. 



3. X states with the identical local purity 

We call a 2-qubit state the X state if the only nonzero 
elements in the density matrix lie along the diagonal or 
skew diagonal. By local unitary transformations, the en- 
tries in the R matrix can always be written as 



x = (0,0, X3), 
T = diag(ii,t2, 



ts). 



(0,0,2/3), 



respectively. The coupling strength k of the AD channel 
is set to be 0.02. We calculate G{p^^{t)) and G{pf^'{t)) 
(or simply Gk{t) and G'^(t)) with fc = I, I, and plot the 
results in Fig. 1. We observe several nontrivial aspects. 

(i) The time evolution of quantum correlation, whether 
for system qubits AB or for the environments A'B', is 
neither monotone nor smooth. There exist critical times 
on which the evolution rate exhibits a sudden change 
in behavior. The sudden change may even happen twice. 
More importantly, the quantum correlation in the system 
state can increase during the time evolution, even if the 
state p"^^ is suffering a dissipative process. 

(ii) By comparing Gi{t) with G[{t), we see that there is 
a symmetry between them. Roughly speaking, the evolu- 
tion of Gi from i = to i — > 00 is similar to the "inverse" 
evolution of G}, namely from t -^ 00 to t = 0. Precisely, 
we have the following relationship: Gi(0) = Gj(-l-oo), 
Gi(ii) = Gi(tl), Giih) = G'M), and Gi(+^) = Gi(0), 
where ti and t^ {i = 1,2) are the critical time for Gi{t) 
and Gj(t), respectively (see Fig.l). In other words, there 



(I) 




(11) 



020 - 
OlS - 
010 - 


WhC) 


\ 








- - 


,, , 



sure, out of a CC state. Our work provides a concrete 
evolution process, in particular, with respect to the state 
Pj^, in which we see a rapid increase of the quantum 
correlation in system qubits. It should be noted that 
the input state p^^ is a separable but not a CC state. 
Meanwhile, we also see that the quantum correlation of 
the environments undergoes a similar evolution: Gj(t) 
increases in the beginning and then changes suddenly to 
decrease. 

It should be mentioned that quantum discord can in- 
crease in the non-Markovian environments [35], whereas 
the AD channel corresponds to a Markovian process. 

Concerning the asymptotical behavior of G{t), we see 
that G{t) oc 7^(i) — e^'"* after the second sudden change 
(if existing). This behavior is qualitatively identical to 
the decay of the skew diagonal elements of the den- 
sity matrix p^^{t), e.g., (p"^'^(i)) . ex e^'^*. It means 
that the asymptotical evolution of quantum correlation is 
closely related the coherence decay (in the case discussed 
above, they are qualitatively identical). Our discussion, 
from the viewpoint of symmetric distance quantifier, pro- 
vides an evidence for the robustness of quantum discord 
[•34] and the claim that "almost all quantum states have 
nonclassical correlations" [42]. 



FIG. 1. Plots of the quantum correlation evolutions G{t) 
(solid line) and G'{t) (dashed line) for states p'j^ (upper 
panel) and p j^ (lower panel) respectively. See text for de- 
tails. 



is a correspondence at the initial, critical and asymptot- 
ical point. A similar situation occurs in the comparison 
of G'ii(i) andG';i(i). 

(iii) The fact that G{0) = G'{+oo) means that the 
initial quantum correlation between the system qubits 
is transferred completely to the asymptotical quantum 
correlation between the ancillary qubits. 

Let us make some remarks about the above observa- 
tions and related topics. 

For a 2-qubit system interacting with two independent 
environments, the dynamics and the sudden change be- 
havior of quantum discord have been analyzed in Ref. 
[32, 33]. Subsequently, these issues are discussed in [39] in 
terms of another quantifier of quantum correlation, one- 
sided HS distance measure [23]. The results we present 
here are referring to the two-sided (i.e., symmetric) HS 
distance measure. One significant phenomenon is that 
quantum correlation G can increase under the influence 
of the separate independent dissipative environments. 

Recently, some authors show that the quantum discord 
can increase under a local amplitude damping channel 
with a CC input state [4(J]. A more general conclusion 
is drawn in [41], that is, any local channel which is non- 
unital and not semi-classical can in principle create quan- 
tum correlations, independently of the considered mea- 



V. CONCLUSION 

In conclusion, we have introduced in this paper a sym- 
metric geometric measure of quantum correlation. By 
performing two-sided von Neumann measurements on bi- 
partite state, we obtain a MICC state. The geometric 
measure is defined as the HS distance between the given 
state and the closest MICC state. For 2-qubit system, we 
simplify the optimization procedure considerably, that 
is, the two-sided optimization is reduced to the one-sided 
one. Hence the numerical computation can be performed 
efficiently. Moreover, the analytical results are available 
for some special class of 2-qubit states. 

Using this quantifier, we have studied the dynamics of 
quantum correlation under the action of AD channel. We 
present the nontrivial aspects which may be exhibited 
during the time evolution: (i) the quantum correlation 
can increase; (ii) the quantum correlation can change 
suddenly, even twice. As for the environments, we see 
that the quantum correlation therein increases from zero 
at the beginning, and then evolves asymptotically to the 
value of the initial quantum correlation in the system 
qubits, after one or two sudden changes. 

The geometric measure developed in this paper pro- 
vides a symmetric viewpoint to study the quantum cor- 
relation. It is also a measurement-oriented measure, since 
the minimal HS distance is referring to the MICC, rather 
general CC, states. Then a question follows: Is there 
more rigorous quantifier, which is referring to the gen- 
eral CC states? This question has been solved in the 
case of one-sided geometric measure [^)N]. However, it 
remains open for the two-sided measure. 
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Appendix A: G for X state with the same local 
purity 

Let \x^\ — IJ/3I = r with < r < 1. We assume that 
all ti is nonzero for the sake of simplicity. Referring to 
(27), we rewrite XM{t) as 



AAf(r 



r^+l'^ + J (r2 - f 2)2 + i^rHlil 



+ r^e^ 



Define the function F{£ ) as 



3- 

(Al) 



(A2) 



In the following, we discuss two cases separately: ti — t^ 
and ti 7^ ^2- 



Case of ii = f2 



When ti=t2, we have f ^ = ^2 _ (^q _ ^2)^2^ g^ ^^^^ 
functions F{i) and Xm{^) depend only on ^§. By taking 
derivative of Am (^3) with respect to £|, and solving the 



equation 



d\M{f4) _ 



for ^3, we have the following results. 



If tits > and r^ G [^3(^1 — ^3), ti{ti — t^)], we have 



*3 



-r^+ti{ti~t3) 

(ti - hr ' 



(A3) 



If iits < and r^ e [-^3(^1 +^3), ^1(^1 +^3)], we have 



*3 



-r^+tijti+ts) 

ih + ts)^ ■ 



(A4) 



Some remarks are needed here. 

(i) The conditions for (A3) and (A4) come from two 
considerations: one is the non- negativity of F{i'^); the 
other is the requirement of f| S [0, 1]. 



(ii) In solving the equation 



dXMiei) 



0, it is assumed 



d{q) 

that F{el) ^ 0. In fact, if F{el) = 0, we have ^3 = 
and r^ = £' 2 . It follows that Xm = r^ ■ We will see below 
that r^ cannot be the maximal value of Aa/(^|). 

(iii) Eqs. (A3) and (A4) require that ti ^ t^ and ii 7^ 
— ^3, respectively. In fact, if ti = ^3, we have ii = ^2 = 
t3 = t and f 2 ^ t2. It follows from (Al) that 



^Min^l^ 



r^ + t^ 



{r^-f^Y+ArH'^il 



+ r^£ 



Inl 



It is a monotone increasing function for £3 € [0, 1]. Then 
A^'*'^ = 2r2-f ^2 when il = 1- This result can be contained 



in (A8). Similarly for the case of ti — —t^. This analysis 
also holds for the derivation process in the next Case. 
Inserting (A3) and (A4) into (Al) respectively, we get 



two candidates for A™/'^: 



(1) _ r2[2ti(^i-t3)~r2] 

(ii - tsr 

(2) _ r2[2ti(fi+t3)-r2] 

M ^ " 



{h+tsr 



(A5) 
(A6) 



We have to take the end points of ^§ into consideration, 
i.e., il = and i^ = I. When i^ — 0, we have Xm = 



-ti\\/2. Let's define 



A 



(3) 
M 



max{r , ti} 



When ^3 = 1, we have 



A 



(4) 
M 



2r2 + ti 



(A7) 



(A8) 



Combining the above analysis, we conclude as follows. 

Ifiits > and r2 g [tsiti—t^), ii(ti—i3)], the maximal 
Am(^) is given by 



\ max 



max{Aiy,2r2 + t2,t2|^ 



(A9) 



If iii3 < and r2 e [-i3(ti + ^3), ti{ti + ^3)], the 
maximal Xm {£ ) is given by 



x^r = ^^^{^^^,^r^+tltl}- 



2. Case of ti / t2 



(AlO) 



Let's prove the following lemma. 
Lemma 2. Ifti 7^ t2, then at least one of ii is zero. 

With Am(^) given by (Al), we introduce Lagrange 
multiplier ^ and take partial derivative of XM{i ) +m(^^ ~ 
1) with respect to £{. Three equations follow. 



2/^2 



tiir 



PI-2\ 



F 

a tlir^-n 



F 



F 



2^ 
2[i 
2r2 



= 0, 

= 0, 

1[i 



0. 



(All) 
(A12) 
(A13) 



Assume that all ii are nonzero. If so, we can delete ii 
in each equation. By noting that all ti are nonzero (as 
assumed at the beginning of this Subsection) , we see from 
(All) and (A12) that the Lagrange multiplier /i must be 
zero. Then Eqs. (All) or (A12) reduce to 1 - {r^ - 
i''^)/\/F = 0. It follows that ^3 = 0, which contradicts 
the assumption. Lemma 2 is proved. 

Subsequently, we will discuss one by one the measure- 
ments allowed by Lemma 2. In each case, we obtain a 
candidate for A™/''. The largest one is what we want. 



a. If £i ^ and £2,^3 7^ 0. It follows that f ^ ^ 
*2 - (*2 - tD^h and that both \m0) and F(/) are the 
functions of ^| only. By the approach similar to that 
presented in Section A 1, we have the following results. 

If t2t3 > and r^ e [t3{t2 - h), ti{t2 - h)], we have 



,(5) _ r^[2t2{t2 - ts) - r^] 



{t2 - hr 



(A14) 



d. If £3 : 

follows that 



1. In this case, 



,(10) 



Aflf — 2r + <3. 



0. It easily 



(A19) 



It is not difficult to see that the above four cases cover 
all allowed measurements. We summary the results ob- 
tained in this Subsection in Table I. 



If <2t3 < and r^ e [-^3(^2 + ^3), ti(i2 + is)], we have 



.^6) _r^[2t2{t2+t3)-r^] 



M/ 



{t2 + hy 



(A15) 



b. If £2=0 and h,h + 0. In this case, £'^ ^t{- 
{t\ — t^£\. The results are very similar to that presented 
in the last paragraph, that is: 

If ii^a > 0, we have 



\(^) 

^M 



r^[2hiti - ts) - r^ 



ih - hy 



(A16) 



If tits < 0, we have 



(8) _ r^[2ti{ti + ts) - r^ 
iti+t^r 



(A17) 



Note that Xj^J and X]^J have the same form as X\j and 
X\,/ (given in (A5) and (A6)) respectively. 

c. If£3 — 0. Here it is not required that both £2 and 
£3 are nonzero. Inserting £3=0 into the expression of 
''^m(^)j i-e., Eq. (Al), we have 



Now let's use an example to show how to calculate 
A™/'^. Given an X state p^^ with |a;3| — ly^] — r, we 
write its R matrix. If we see that ^2^3 > and ^1^3 > 0, 
then we face the following possibilities: 



(i) If r^ G Intervali and r^ e Interval3, then A™/'^ 
raajc\Ajyj ,Aj^ ,A^j ,Aj^j /• 

(ii) If r^ e Intervali and r^ <^ Interval3, then A' 
max-[A^j , Aj^j , Aj^j /• 

TABLE I. Summary of the results in Section IIIB3 



max 

M 



£ 


Conditions 


Am 




£1 = 0, and 
4,4/0 


^2*3 > 0, r^ G Intervali "" 
t2tz < 0, r^ G Interval2 


X(6) 


(A14) 
(A15) 


£2 = 0, and 

4,4/0 


fits > 0, r^ G Intervals 
fits < 0, r^ G IntervaU 


X(8) 


(A16) 
(A17) 


4 = 




X(9) 


(A18) 


4 = 1 




x(i()) 


, (A19) 



'^ Intervali = [tz(t2 - t3),t2(t2 - ta)], 
Intervals = [-t3(t2 + ts), t2{t2 + ts)], 
Intervals = [tz(ti - tz),t-i{ti — ts)], 
IntervaU = [-tsC*! + ts), ti(ti + ts)]- 



Am(/) = \v^ + i4el + 4^1) - \r^ - itl£l + tl£l) 



The maximal value of the above expression is given by 



(iii) If r^ ^ Intervali and r^ G Intcrval3, then A 

,(7) x(9) >(!( 



„-,„„r\(7) \(9) ^(W)l 

max-[A^ , Aj^^ , Ajy^ /• 



max 

M 



(iv) If r^ ^ Intervali and r^ ^ Intervals, then A 
max{AJ{^\A^^°^}. 



max 

M 



Ai^^=max{r2,t2,t2}. 



(A18) 



Then quantum discord G{p) is obtained by inserting A 
into (29). 



max 

M 
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